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We have mapped in the 2.7 mm continuum and 12 CO with the PdBI the IR-dark "tail" that crosses the IC 1396N globule from south 
to north, and is the most extincted part of this cloud. These observations have allowed us to distinguish all possible associations of 
molecular hydrogen emission features by revealing two well-collimated low-mass protostellar outflows at the northern part of the 
globule. The ouflows are located almost in the plane of the sky and are colliding with each other towards the position of a strong 
2.12 yum H 2 line emission feature. 
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1. Introduction 

IC 1396N is a bright-rimmed cloud (BRC38; Sugitani et 



al. U991l ). located at a distance of 750 pc (Matthews [1979l l. which 
has been extensively studied by our group (Nisini et al. 2001 1 
C odella et al. [2001 ; Beltran et al. [20Tj2l |20D4l [20091 Fuente et 
al. 2009). The cloud shows a cometary structure elongated in the 
south-north direction, with the globule head located at the south, 
facing the 06.5 exciting star HD 206267, and the tail pointing to 
the north. Towards the globule head, one finds the intermediate- 
mass protostar IRAS 21391+5802 (BIMA 2), which is asso- 
ciated with multiple millimeter compact sources and powers a 
molecular outflow (Codella et al. 2001 ; Beltran et al. 2002 ; Neri 



et al. 2007; Fuente et al. 2009 ). Near-infrared (NIR) images have 
revealed a number of small-scale 2.12 fim H2 line emission fea- 
tures and Herbig-Haro (HH) flows spread all over the globule 
(Nisini et al. |2001; Sugitani et al. 125021 Reipurth et al. 120031 
Caratti o Garatti 2006; Beltran et al. |2009l , as well as an IR-dark 
"tail" that crosses the cloud from south to north and is clearly 
pinpointed in the JHK' composite image (see Fig. 1 of Beltran 
et al. |2009l and Fig.[Qi). This IR-dark "tail" is well-traced in both 
CO and high-density tracers such as CS and H 13 CO + (Codella et 
al. l200Tl Sugitani et al. l2002b . 

Deep sub-arcsecond 2.12 fim H2 1-0 S(l) line emission 
observations carried out with the TNG telescope (Beltran et 
al. 12009 ) resolved the emission features spread throughout the 
region into several chains of knots. In a few cases, these strands 
of H2 knots have a jet-like morphology, which suggests that 
they could trace different flows (Nisini et al. 120011 Beltran et 
al. 120091 1. Towards the northern part of the globule, Nisini et 
al. {2001) proposed that some of these knots could be associ- 
ated with the CO molecular outflow discovered by Codella et 
al. (120011) . and Beltran et al. d20091 > suggested that strand of 



knots located at both sides of the IR-dark "tail" could be as- 
sociated with each other and be part of long flows with lengths 
ranging from 0.18 to more than 0.4 pc, and that the powering 
sources were embedded in the "tail". 

To test the hypothesis that some of the molecular hydrogen 
emission features could be associated with long flows cross- 
ing the northern part of IC 1396N, we mapped with the IRAM 
Plateau de Bure Interferometer (PdBI) the whole IR-dark "tail" 
in 12 CO, which is the most sensitive outflow tracer and should 
help us to trace the cold-warm flow component that is over- 
looked by the hotter component tracer H2. The "tail" was also 
mapped in the continuum at 2.7 mm to search for the possible 
outflow powering sources that according to Beltran et al. (2009 ) 
should be hidden there. In this Letter, we present the first results 
of this study, which confirm the scenarios proposed by Beltran 
et al. (2009 ), including the collision of two of these flows. 

2. Observations 

We carried out seven-point mosaic observations of the IR-dark 
"tail" in IC 1396N in the 2.7 mm continuum and 12 CO (1-0) line 
emission with the PdBI on 2010 July 24 and 27, and November 
4, in the D and C configurations. The inner hole in the (u, 
v)-plane has a radius of 7 Li. Line data were smoothed to a 
spectral resolution of 1 kms~'. The maps were created with 
natural weighting, attaining a resolution of 2'.' 9 x 2'.'8, but to 
highlight the continuum and line emission, the final maps were 
reconstructed with a circular synthesized beam of 4" x 4". 
The lcr noise is 0.25 mJybeanT 1 for the continuum map and 
10 mJybeatrT 1 channeL 1 for the CO line map. 

3. Analysis and discussion 
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The continuum observations of the northern part of the IR-dark 
"tail" detected two embedded cores, which we call sources I and 
C, powering well-collimated outflows (with collimation factors 
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Fig. 1. (a) H-band image of IC 1396N taken with NICS at TNG. The white crosses show the positions of the 3.1 mm sources, BIMA 
1, 2, and 3 from Beltran et al. (12002V the Class I sources # 252 and # 331 (bina ry, A+B) (Beltran et al.[2009), the tw o embe dded 
sources I and C, which are powering the outflows in panel b (Codella et al. 1200 It this paper), and core II (Sugitani et al. 2002). The 
dotted and dashed boxes indicate the region mapped in panels b and c, respectively, (b) H2 (2.12 //m) image (continuum-subtracted) 
in grey-scale (Beltran et al. 2009 1, 2.7 mm continuum emission in green contours and 12 CO (1-0) emission averaged over the 
[+8, +38] km s~ velocity interval in red and blue contours towards the northern part of the IC 1396N cloud. Contour levels are 3, 
6, 9, 12, 15, and 24 times lcr, where lcr is 0.25 mJybeairT 1 for the continuum, and 4 mJy beam 1 for the CO emission. The PdBI 
synthesized beam is shown in the lower right corner, (c) H2 (continuum-subtracted) (top), //-band (middle), and /-band (bottom) 
image towards the chain of knots C (Beltran et al. 120091 1 . 



>5 towards the central part of the outflows; Fig. [TJ)) and sepa- 
rated by -0.1 pc. The source C was first detected at 1.25 mm 
with the IRAM 30-m telescope by Codella et al. ( 120011) . who 
suggested that it was the powering source of the northern bipo- 
lar outflow mapped in 12 CO (2-1) with an angular resolution of 
10". However, this is the first time that the source is mapped 
and at high angular resolution. The existence of source I and the 
bipolar outflow that is driving had been suggested by Beltran et 
al. ((2009 ) (see next section), and our observations here confirm 
this idea. This new source could be associated with the H 13 CO + 
peak labeled as core I by Sugitani et al. (2002). The mass of the 
sources was estimated assuming a dust opacity of 1.15 cm 2 g~' 
at 1 15 GHz (Beckwith et al. 1990) and a gas-to-dust ratio of 100. 
Since the dust temperature, r^st, of the sources is unknown, we 
estimated the masses for 10 K and 20 K, which are the typi- 
cal temperatures of the gaseous envelopes of low-mass proto- 
stars. The mass of source I is -0.5 M (for 7d us t=10 K) and 



-0.2 M (for r dust =20 K), and that of source C is -0.3 M (for 
7du.st=10 K) and -0.1 M (for r dust =20 K). 

3.2. The colliding molecular outflows 

The elongated structure and the very well-defined bipolar mor- 
phology of the two molecular outflows mapped in 12 CO (1-0) 
(Fig. QJ>), suggest that they are close to the plane of the sky. 
Both outflows are associated with 2.12 /mi H2 emission. The 
blueshifted lobes of the two bipolar outflows strongly overlap 
towards the position of the H2 strand of knots C. The northern 
outflow, named N, is being powered by source C and is ori- 
ented in the northeast-southwest direction. Nisini et al. (2001) 
and Beltran et al. (2009) suggested that it could be associated 
with the strand of H2 knots D and C. As seen in Fig. [TJi, the 
CO emission indeed complements the H2 emission where the 
excitation conditions and the opacity of the material are differ- 
ent. On the other hand, the southern outflow, named S, which is 
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Table 1. Properties of the molecular outflows" '' 
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(a) Estimated values. The outflow parameters should be corrected if the inclination angle with respect to the plane of the sky, i, is known. 

(b) Blueshifted emission averaged over the velocity range [-38, -8] km s -1 and redshifted one over [+8, +38] km s" 1 . 

(c) Momenta and kinetic energies are calculated relative to the cloud velocity, Vlsr=0 km s . 



driven by source I and oriented in the southeast-northwest direc- 
tion, had never been mapped before. Beltran et al. (2009 1 had 
proposed that the chains of H2 knots F, T, and C could be as- 
sociated with and be part of the same flow, which would have a 
length of ~0.4 pc, as the CO observations have confirmed. These 
authors had also predicted that the driving source, in this case 
source I, should be embedded in the dense gas mapped in CS 
by Codel la et al. (T200TT) and in H 13 CO + (core I) by Sugitani et 
al. (120021 1. According to the correlations between driving source 
and outflow properties of Cabrit and Bertout (1992 1, the bolo- 
metric luminosity of the powering sources should be < 10 L Q . 
The parameters of the outflows, which are estimated assuming 
that the CO emission in the wings is in local thermal equilibrium 
and optically thin, are given in Table 1 . The outflow parameters 
and the masses of the driving sources are consistent with those 
of low-mass protostars (e.g., Arce et al. 120071 and references 
therein). Since the CO spectra exhibit prominent self-absorption 
at the cloud velocity and it is impossible to derive the excitation 
temperature, r ex , from the brightness temperature, the outflow 
parameters were estimated assuming r ex =10 K, which was the 
value adopted by Codella et al. ([2001 ) for the northern outflow. 
We used an Ithl/ICO] abundance ratio of 10 4 (e.g. Scoville et 
al. 119861 1. To calculate the parameters of the blueshifted lobes, 
we assumed that half of the CO emission seen towards the fea- 
ture C is associated with outflow N and half with outflow S. That 
the parameters of the redshifted lobes are similar to those of the 
blueshifted lobes for both outflows suggests that this assumption 
is correct. The parameters need to be corrected by the factors 
indicated in Table 1 if the inclination angle of the flow with re- 
spect to the plane of the sky, i, is known. On the basis of the 
morphology of the outflows, i < 10°. 

In the scenario proposed by Beltran et al. (2009), the out- 
flows N and S would collide towards the position of the chain of 
knots C and this would explain the strong H2 emission observed 
towards this feature. As previously mentioned, the blueshifted 
lobes of the CO bipolar outflows do indeed overlap towards the 
position of the strand of knots C (Fig.QJ). However, the collision 
of outflows is not the only possible explanation of this morphol- 
ogy. The overlapping of the blueshifted emission might also in- 
dicate that the outflows are colliding with the same dense clump, 
or are merely a projection effect. In the latter scenario, the out- 
flows would be in different planes of the sky and would be shock- 
ing with different dense clumps aligned along the line of sight. 
However, that the outflows are almost in the same plane (that 
of the sky), that the powering sources are at the same systemic 
velocity, Vlsr~ km s , and that some observational facts fa- 
vor a collision (see next), either with the same clump or between 
outflows, rules out a chance superposition of the outflows. 



Figure [2] shows the position-velocity (PV) plots along the 
major axis of the outflows. The cuts were made along PA=- 102° 
for outflow N and PA=-24° for outflow S, and all the plots were 
obtained after averaging the emission along the direction per- 
pendicular to the cut for the purpose of increasing the S/N of 
the plots. The clumpy morphology of the outflows (Figs.flJ) and 
13 could be the result of episodic mass-loss events. As seen in 
Fig. |2l the velocity of the CO emission increases rapidly with 
increasing distance from the powering sources. This is clearly 
visible for the southern outflow. The material accelerates until 
it reaches the position of the H2 knots: the strand of knots D 
and F, located at the end of the redshifted lobes of the N and S 
outflows, respectively, and the group of knots C, at the end of 
the blueshifted lobes. At the position of the H2 knots, there is 
an increase in the extreme emission velocity (V>40 km s _1 ) for 
the lowest emission contour level, while the bulk of the emis- 
sion abruptly decelerates and has a much lower velocity ( V^IO- 
20 km s~'). This is less evident for the knots D, which reach 
velocities of only ~20 km s -1 , although this could be due to the 
lower sensitivity at the edge of the map to higher velocities. For 
the chain of knots C, the bulk of the emission clearly delineates 
an inverse triangle shape for both outflows, with the velocity of 
the material slowly increasing from the position of knot C5 (the 
first knot of the main H2 emission feature encountered by the 
blue lobes) to the position where the material at the lowest emis- 
sion level reaches the maximum velocity, and then decelerating 
up to the position of knot C16 (the last knot of the main H2 emis- 
sion feature encountered by the blue lobes). The morphology of 
the CO emission could be explained in terms of outflowing ma- 
terial that is impacting, either with a dense clump or with another 
outflow, at the position of the H2 knots, giving rise to the NIR 
emission. 

The H2 emission of the group of knots C is the strongest of 
the entire northern part of the globule. The integrated line emis- 
sion of this feature, summing the emission of all the individ- 
ual knots, is 485.4xl0~ 15 mWirT 2 (Beltran et al. |2009l . This 
value is much higher than the integrated molecular hydrogen 
emission of the strand of knots D (52.8xl0~ 15 mWirT 2 ) and F 
(92.9xl0~ 15 mWirT 2 ). This clear asymmetry at both ends of the 
outflows, in terms of both morphology and intensity of the fea- 
ture, led us to think that the shock that is taking place at the posi- 
tion of the knots C is exceptional. The knots D, F, and C are also 
visible at 4.5 /mi (see Fig. 5 of Beltran et al. 2009), in which case, 
the feature C is also the strongest. We investigated whether this 
could be due to extinction. Assuming that knots D and F have 
the same integrated emission as knots C but are more extincted, 
to reproduce the same integrated fluxes, the visual extinction, 
Ay, should be ~21 mag higher in D than in C, and ~16 mag 
higher in F than in C. From the color-magnitude of Beltran et 
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al. (2009 '), and from estimates based on the CO and CS emission 
data of Codella et al. (12001 1. the maximum extinction through 
the northern end of the globule "tail" is Ay < 20 mag. Although 
we cannot exclude that the extinction is severely affecting the 
intensities of the knots D and F, it seems unlikely that the differ- 
ence in extinction between molecular hydrogen features located 
outside the IR-dark "tail", which is the most extincted part of the 
globule, are so large. Spectra taken towards the position of knots 
F and C (Massi et al. in preparation), show a richness of other 
H2 transitions in J-, H-, and /T-bands. The emission is much 
stronger in C than in F, which would explain why the strand of 
knots C is clearly visible in the H and /-band images (Fig.QJ). 

The CO and H2 emission clearly indicates that there is a 
shock at the position of the strand of knots C. To distinguish 
between the two possible scenarios, that is, that the outflows are 
colliding with each other or are impacting with the same dense 
clump, we analyzed the emission at the cloud velocity traced by 
12 CO in the [-8, +8] km s velocity range (the cloud Vlsr is 
~0 km s _1 ). The average CO emission in this velocity interval 
does not show any compact clump at that position. The nega- 
tive contours visible in the PV plots at the position of the knots 
C at systemic velocities would indicate that the interferometer 
has filtered out only extended emission. A compact dense clump 
is similarly invisible in the channel maps of the high-density 
tracer CS (3-2) (Codella et al. 1200 11 1. Therefore, the most plau- 
sible explanation for the NIR and CO emission is that the out- 
flows are located in the same plane, as already suggested by their 
well-defined bipolar morphology, and are indeed colliding. This 
scenario is particularly interesting because the collision would 
take place in a very pristine environment, where sources C and 
I would be the only embedded protostars. As seen in Fig. QJ, 
most of the protostars are located in the southern part of the IR- 
dark "tail". As far as we know, possible outflow interactions have 
previously been reported in the crowded LI 448 region, where 
multiple protostars and molecular outflows have been detected 
(Barsony et al. 119981 Kwon et al. 120061) . However, the most 
likely explanation of those interactions are not outflow collisions 
but outflow-dense core shocks (Barsony et al.[T998). 

Assuming that both outflows are randomly oriented, the 
probability of such an event is low (~5%), if one considers the 
case of two sources separated by a short distance (0.2 pc), power- 
ing outflows with 0.01 pc of radius, in which one of the outflows 
is perpendicular to the line joining the two sources. The prob- 
ability is higher (~ 10-20%) if one considers the case in which 
the direction of both outflows is allowed to vary. We note that 
the probability will be even higher if the separation between the 
sources is smaller or the opening angle of the outflow higher. 
Therefore, that, to our knowledge, clear cases of outflow colli- 
sions have not been found is not because of the rather low proba- 
bility of such an event but to the practical difficulty in observing 
two very close sources in a similar evolutionary state located at a 
distance that allows us to resolve the emission of their respective 
outflows. 
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Fig. 2. 12 CO (1-0) position-velocity plots along the major axis 
of the outflows, with PA=-102° for outflow N, and PA=-24° 
for outflow S. The offsets are measured from the position 
of source C, ar(J2000)=21 h 40 h 39 s .04, 5(J2000)=+58°18'29'.'8, 
for outflow N and of source I, a(J2000)=21 h 40 h 40 s .66, 
<5(J2000)=+58°18'06'/0, for outflow S, and are positive towards 
the northeast and southeast, respectively. Contour levels are -12, 
-9, -3, 3, 6, 9, 12, 27, 45, 90, and 145 times lcr, where cr is 
10 mJybeairT 1 . The horizontal line indicates the Vlsr- The ver- 
tical dotted line indicates the position of the powering sources, 
and the vertical dashed lines the position of the H2 knots (Beltran 
etal. 120091 . 

Additional high-angular resolution observations of shock trac- 
ers, such as SiO, would be needed to confirm this collisional 
scenario. Observing lower SiO velocities over the blueshifted 
interacting region than for the redshifted emission of the out- 
flow would further support this scenario. The collision scenario 
would imply, if confirmed, that the star formation in this globule 
is aligned along a privileged direction. 



4. Conclusions 

We have presented evidence for the collision of the outflows N 
and S observed in the northern part of the bright-rimmed cloud 
IC 1396N based on the overlap of the CO blueshifted lobes, the 
well-defined bipolarity of the outflow, which indicates that both 
outflows are located close to the plane of the sky, the strong 
H2 emission towards the position where the two blueshifted 
lobes overlap (strand of knots C), and the non-detection of a 
compact dense clump at the position of the putative shock. 
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